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Abstract 



I We consider elastic reflection and transmission of electrons by a disordered 
system characterized by a 2N x 2N scattering matrix S. Expressing S in 

■ terms of the radial parameters and of the fom' NxN unitary matrices used 
. for the standard transfer matrix parametrization, we calculate their prob- 
, ability distributions for the circular orthogonal (COE) and unitary (CUE) 
I Dyson ensembles. In this parametrization, we explicitely compare the COE- 

■r^lj- ' CUE distributions with those suitable for quasi-ld conductors and insulators. 

■ Then, returning to the usual eigenvalue-eigenvector parametrization of S, we 
"j^ ■ study the distributions of the scattering phase shifts. For a quasi-ld metallic 
^ . system, microscopic simulations show that the phase shift density and corre- 

I lation functions are close to those of the circular ensembles. When quasi-ld 

I longitudinal localization breaks S into two uncorrelated reflection matrices, 

^ ■ the phase shift form factor b{k) exhibits a crossover from a behavior charac- 
teristic of two uncoupled COE-CUE (small k) to a single COE~CUE behavior 



(large k). Outside quasi-one dimension, we flnd that the phase shift density 
is no longer uniform and S remains nonzero after disorder averaging. We use 
I perturbation theory to calculate the deviations to the isotropic Dyson distri- 

butions. When the electron dynamics is no longer zero dimensional in the 
transverse directions, small-A: corrections to the COE-CUE behavior of b{k) 
appear, which are reminiscent of the dimensionality dependent non universal 
regime of energy level statistics. Using a known relation between the scatter- 
ing phase shifts and the system energy levels, we analyse those corrections to 
the universal random matrix behavior of S which result from d-dimensional 
diffusion on short time scales. 
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I. INTRODUCTION 



The discovery of universal conductance fluctuations (UCF characterizing the sensi- 
tivity of the conductance of small metallic samples to a change of the Fermi energy, magnetic 
field or impurity configuration) have generated a sustained interest in quantum mesoscopic 
physics since the mid eighties. Mesoscopic systems have a size of the order of the electron 
phase-coherence length L^, i.e. a scale intermediate between single atoms (microscopic) 
and bulk solids (macroscopic). First, it is in terms of quantum interference effects between 
different multiple scattering paths that UCF has been understood |@J^. The universality of 
the phenomenon (reproducible fluctuations of order e'^/h, independently of the mean con- 
ductance) quickly lead physicists 0,^ to understand UCF as a signature of a more general 
universality resulting from the eigenvalue correlations of random matrices. 

The standard random matrix ensembles were introduced in the context of Nuclear Physics 
and later found to also describe the statistical properties of quantum systems whose 



classical analogs are chaotic [p!0| , pA| . Contrary to a microscopic approach where the system 
hamiltonian H, scattering matrix 5* or transfer matrix M result from a more or less arbitrary 
distribution of the substrate potential, random matrix theory (RMT) assumes for H, S or 
M statistical ensembles resulting from an hypothesis of maximum randomness, given the 
system symmetries (time reversal symmetry, spin rotation symmetry, current conservation) 
plus a few additionnal constraints. 

A direct relationship between conductance fluctuations and RMT comes from the Thou- 
less expression of the conductance 

G = j N{E^) , (1.1) 

given by the number of one-electron Q levels N{Ec) which lie within an energy band of width 
Ec = Dh/L^ centered around the Fermi energy Ef. The Thouless energy Ec is the inverse 
characteristic time for an electron to diffuse through the sample (of size L). The electron 
diffusion coefficient is D = Vfl/d, where f/ is the Fermi velocity, / the elastic mean free path, 
and d the spatial dimension of the sample. Fluctuations in the number of levels within E^ are 
therefore related to conductance fluctuations. The analysis of fluctuations in the spectrum 
of non interacting electrons in metallic particles was initiated by Gorkov and Eliashberg 
. The relevance of RMT was proved by Efetov |13[ and complemented by Altshuler and 



Shklovskii who noticed non RMT behavior for energy separations larger than Ec- Using 
a diagramatic perturbation theory for the density-density correlation function in the weak 
disorder limit kl ^ 1 {k is the Fermi wave vector), they showed that the correlation between 
energy levels is correctly described by the universal RMT laws for energy scales smaller than 
Ec, but are weaker (and dimensionality dependent) beyond -Ec-Q 



^Throughout this work we will treat spinless electrons and therefore we will not write spin- 
degeneracy factors. 

^This limit of the universal RMT laws is specifically derived ignoring electron-electron interaction. 
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Pionieered by Imry an alternative approach relating UCF and random matrix the- 
ories is based on the Landauer formula, which gives the conductance in terms of the total 
transmission coefficient of the sample (at the Fermi energy), considered as a single, complex 
elastic scatterer. For a two-probe measurement where the disordered sample is attached 
between two perfect reservoirs with an infinitesimal difference in their electrochemical po- 
tentials, the conductance (measured in units of e^/h) is 

9 = Tr[ttt] . (1.2) 

The transmission matrix t can be expressed in terms of the transfer matrix M, for which a 
standard random matrix theory ("global approach") has been developed [l^Jl^. A set of 



real positive parameters describing the radial part of the 2Nx2N transfer matrix M (precise 
definitions are given in the next section) are the relevant "eigenvalues" in this approach, and 
their probability distribution is 

P({AJ)=exp(-/37^({AJ)) , (1.3) 



N N 

^({Aa}) = - E log |Aa - Afel + E V{\a) . (1.4) 
a<b a=l 

This is the usual RMT Coulomb gas analogy with logarithmic pairwise interaction, a system 
dependent confining potential ^(A), and an inverse temperature /? = 1, 2, 4 depending on the 
system symmetries. This distribution corresponds to the most random statistical ensemble 
for M given the average density p(A) of radial parameters, which controls the average con- 
ductance (g). In this maximum entropy ensemble, V{X) and p(A) are related by an integral 
relation in the large A^-limit 

V{\) + C= dA'p(A')log|A-A'| + ^logp(A) . (1.5) 
Jo 2p 

C is a constant, and to leading order in this mean field equation expresses the equilibrium 
of a charge at A resulting from its interaction with the remaining charges and the confining 
potential. 

In this work we examine the statistical properties of another matrix related to quantum 
transport in disordered systems: the scattering matrix 5*. 

First, the relation between M and S is straightforward, since M can be expressed in terms 
of the refiection and transmission submatrices which define S. This allowes us to show that 
the A-statistics characterizing the Dyson ensembles coincides with a "global approach" for 
M, given a particular confining potential V{\). Using this equivalence, the detailed proof 
of it being given in this work, the quantum transport properties associated with the circular 



ensembles have been obtained IjT^, including the weak-localization effects and the quantum 
fluctuations of various linear statistics of the A-parameters (conductance, shot-noise power, 
conductance of a normal-superconducting microbridge, etc) . This equivalence has been also 
derived by Baranger and Mello []TB|, and has been conflrmed in their numerical simulations 
of chaotic billiards. 
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Second, since 5" is unitary, its 2N eigenvalues {exp {iOm)} are just given by the 2N phase 
shifts {Om}- In the Dyson circular ensembles, where all matrices S with a given symmetry are 
equally probable, the phase shift statistics follow the same universal level correlations than 
in the gaussian ensembles of corresponding symmetry. The applicability of Dyson circular 
ensembles for chaotic scattering has been established by the pioneering work of Bliimel and 
Smilansky ||19|, as far as the two-level form factor for the phase shifts is concerned. Further 
evidence comes from numerical studies of 2 x 2 matrices describing the scattering through 



chaotic cavities 20 



But a ballistic chaotic cavity, with a conductance of order N/2 and where the electronic 
motion is essentially zero-dimensional after a (short) time of flight, differs from a disordered 
conductor or insulator. The introduction of bulk disorder in the dot reduces the conductance 
to smaller values of order Nl/L {L >> I , metallic regime) and enhances the time interval 
in which the electron motion depends on system dimensionality. The diffusive motion of the 
carriers after a (short) characteristic time limits (rather than improves) the validity of the 
standard RMT distributions. The subject of this study consists precisely in understanding 
this apparent paradox: the fact that the introduction in a cavity of bulk disorder drives 
the statistics of S away from the circular ensembles. For this purpose, we extensively study 
disordered systems where non-interacting electrons are elastically scatterered by microscopic 
impurities contained in a rectangular dot of various aspect ratios. We derive useful relations 
between different parametrizations of H, M and S, involving the energy levels {Ei}, the 
radial parameters {\a} and the phase shifts {6m}- 

The paper is organized as follows. After presenting the basic definitions of our model 
and the relationship between S and M, we derive (Sec. [llip the metric of 5* in the polar 
decomposition, and therefore the probability distribution of the radial parameters assuming 
a Dyson distribution for the scattering matrix in the time-symmetric case (the case without 
time-reversal symmetry is worked out in Appendix 0). In Sec. |3 we use this parametriza- 
tion to discuss the differences and similarities between the circular ensemble distributions 
and those suitable for quasi-one dimensional disordered systems. In Sec. we start the 
statistical analysis of the phase shifts of S for long conductors with weak disorder. In Sec. 
|VI| we study the scattering phase shifts in quasi- Id insulators and show that their corre- 
lations are well described by those of two uncoupled circular ensembles. In Sec. [VI I| we 
calculate by diagrammatic perturbation theory the mean values of the S'-matrix, which are 
needed to understand the phase shift anisotropy found outside the quasi-lc? limit. In Sec. 



Vlll| we present numerical results for geometries other than quasi-one dimensional strips (i.e. 



squares and thin slabs) and show that after the spectrum is unfolded there is a rough agree- 
ment with the circular ensembles, though noticable discrepancies are visible in the number 
variance. These deviations are analyzed in Sec. p(| , using the known deviations of metallic 
spectra from the R.M.T. behaviors and exploiting a relationship between scattering phase 
shifts and the energy levels (Appendix 0). Appendix ^ gives a summary of the methods and 
results of the numerical simulations, and in Appendix 0, we use a semiclassical approach to 
complement our understanding of the results, including the meaning of the Wigner time. 
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II. SCATTERING, TRANSFER, TRANSPORT AND A-PARAMETERS 



Using only the system symmetries (current conservation, time reversal symmetry, spin 
rotation symmetry) one can show that both S and M can be expressed in terms of the N 
radial parameters {Aa} of M and 4 (2 in the presence of time reversal symmetry) N xN 
auxiliary unitary matrices. In this parametrization, we shall show the similarities and the 
differences between the distributions implied by the circular ensembles and those describ- 
ing in the weak scattering limit long quasi-ld disordered conductors and insulators. To 
make concrete our explanations, we introduce in what follows the essential elements of the 
particular microscopic model on which we test the validity of Dyson circular ensembles. 

We consider an infinite strip composed of two semi-infinite perfectly conducting leads 
of width Ly connected by a disordered part of same width and of longitudinal length 
(Fig. 1). Assuming non- interacting electrons and hard- wall boundary conditions for the 
transverse part of the wavefunction, the scattering states in the leads at the Fermi energy 
Ef — Ti^k'^ /2m satisfy the condition k"^ — (nn/LyY + k"^, where k is the Fermi wavevector, 
nn/Ly the quantized transverse wavevector and kn the longitudinal momentum. The various 
transverse momenta labeled by the index n (n = 1, . . . , A^) which satisfy this relationship 
with k^ > define the propagating channels of the leads. Since each channel can carry 
two waves travelling in opposite directions, asymptotically far from the disordered region the 
wave function can be specified by a 2A"-component vector on the two sides of the disordered 
part. In each lead the first A^ components are the amplitudes of the waves propagating to 
the right and the remaining A^ components are the amplitudes of the waves travelling to the 
left. 

N 1 

^lix, y) = ^-— (A^e'^-^ + 5,e-'=--) </.„(y) , (2.1a) 

n=l Kn 



N 

^>u{x, y) = ^ ^ (C7,e^'="(^-^^) + D^e-'^-^^-"^^^) 4>n{y) • (2.1b) 

n=l kn 



The transverse wavefunctions are 0„(|/) = \J2/Ly sm{7rny/Ly). The normalization is chosen 
in order to have a unit incoming fiux in each channel. The scattering matrix S relates the 
incoming fiux to the outgoing fiux 

S is a 2N X 2 A" matrix of the form 

(2.3) 

The refiection (transmission) matrix r (t) is an N xN matrix whose elements r^a (tba) denote 
the refected (transmitted) amplitude in channel b when there is a unit fiux incident from the 
left in channel a. The amplitudes r' and t' have similar meannings, except that the incident 
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flux comes from tlie riglit. f\ Current conservation implies that S is unitary. Note that, with 
the convention we have taken, S for a perfect (non-disordered) sample at zero magnetic field 
is not the identity matrix but is characterized by transmission submatrices which contain 
pure phases Ua = t'^a = ^a^exp {ihL^). 

The 2N X 2N transfer matrix relates the flux amplitudes on the left-hand side of the 
disorder part with those on the right: 

\ f A \ 

(2.4) 

Just as for S*, we can write M in terms of four NxN blocks 

M = f "^A . (2.5) 

The reflection and transmission matrices of S can be expressed in terms of the block 
matrices rrii of M. Introducing the polar representation pl|j2^ of M, we have: 

r = -m^^m-i = -u^^'^Tlu^^^ , (2.6a) 
t = {m\Y^ = u'^^^Tu^^^ , (2.6b) 
r' = m2mY = u^^^Uu^'^'^ , (2.6c) 

t' = mr^ = u^^^Tu^^^ , (2.6d) 

where m'-'^ (/ = 1, . . . ,4) are arbitrary NxN unitary matrices. TZ and T are real diagonal 
NxN matrices whose non-zero elements (labeled by only one index) are the square roots of 
the reflection and transmission eigenvalues which can be expressed as a function of the real 
positive diagonal elements Aq (a = 1, . . . , N) of the NxN diagonal matrix A characterizing 
the radial part of M: 



In this A-parametrization, M and S can then be written as 



n.-[^) . (2.7a) 



1/2 

= ( ) ■ (2-7b) 



^Throughout this work we shall frequently represent 2N x 2N matrices in terms of their NxN 
blocks. We reserve capital letters for 2Nx2N matrices, while calligraphic and low-case letters are 
used for NxN matrices. 
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M - [ „ ^(2)i M/2 n , \M/2 n „,f3)t ' (2-8) 



Since tt''" = UiT'^Ui, the dimensionless conductance can be expressed as 



TV 

9 



a=l J- + 



The conductance is therefore a hnear statistics of the radial parameters {Xa} of M. Note that 
such a simple relationship does not exists between g and the scattering phase-shifts {Oa}, 
since the relation between the eigenvalues of 5* and those of tt^ depends on the eigenvectors 
of 5* and on the w-matrices. 

In the absence of a magnetic field there is time reversal symmetry, the S-matrix is 
symmetric (5* = S^) and the polar decomposition has only two arbitrary unitary matrices 
since 

u^'^ = , (2.11a) 

n(^) =n(2)T . (2.11b) 

In this case the number of independent parameters is to 2N'^ + N. We have A^^ parameters 
for each of the two N x N unitary matrices and N for the diagonal matrix A. Without 
time-reversal symmetry (unitary case with spin degeneracy) the number of independent 
parameters of S (and M) is 4iV^ (the N extra parameters of the polar decompositions ( p.9|) 
and (|2.8|) are due to the fact that they are not unique |l^). In the symplectic case occuring 
when there is a strong spin-orbit scattering in the disordered part and no applied magnetic 
field, the spin degenracy is removed and each matrix element becomes a 2 x 2 quaternion 
matrix, which doubles the size of M and S, but u^^^ and u^^^ are also given [22] in terms of 
u^^^ and u^^^ and the A have a twofold degeneracy (Kramers degeneracy). 



III. INVARIANT MEASURE OF S IN THE POLAR DECOMPOSITION 

In this section we calculate the invariant measure fi{dS) of S in terms of the radial pa- 
rameters {Aa} and the matrices u^''\ We present here the time-symmetric case (/3 = 1), while 
the unitary case (/? = 2) relevant when a magnetic field is applied is considered in Appendix 
0. Our calculations have recently been extended by K. Frahm |^ for the symplectic case 
{(3 = 4). In the orthogonal case 5* is unitary symmetric and can be decomposed as: 

S = W^E W = U^T U = Y^ Y . (3.1) 

The first equality simply means the diagonalization of S and introduces the 2N phase shifts 
of S through the diagonal elements exp {i6m) of E and a 2N x 2N orthogonal matrix W 
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containing the eigenvectors of S. The second equahty results of the polar representation of 
5* where the real matrix F and the block-diagonal unitary matrix U are given by Eq. ( p.9|) 
with the conditions ( |2.11| ). The last decomposition holds for any unitary symmetric matrix 
and introduces a 2N x 2N unitary matrix Y which is not unique, but specified up to an 
orthogonal transformation. Following Dyson the measure of a neighbourhood dS of 
S is given in terms of the infinitesimal variations dMij of a the matrix elements of a real 
symmetric matrix dM defined by: 

dS = {idM) Y , (3.2) 



2N 

^{dS) = X{dMij . (3.3) 

This definition is independent of the particular choice of the unitary matrix Y and we use this 
freedom of choice to take a convenient Y for expressing dM in the A-parametrization. To this 
end, we note that F is real symmetric and unitary, with eigenvalues ±1 and diagonalizable 
by an orthogonal transformation O: 

F = D O , (3.4) 



--(o • 0-(l %). (3.5) 

The NxN blocks of O are diagonal matrices given by 

= ^l-TZa , 2" = 71 + . (3.6) 

Writing the diagonal matrix D as F^, with 
one can write 5" as Y^Y, with 

*' = ^o£/=(fg"^„ J;;,,). (3.8) 

Since Y is unitary, its infinitesimal variations can be expressed as 

dY = 5YY , (3.9) 

where the matrix 5Y is antihermitic. Analogously, for the bock components of U we can 
write 

du^^^ = , 6u^^^ = da^^'^ + i ds^^'^ , / = 1, 2. (3.10) 
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da^''^ (ds^')) are real antisymmetric (symmetric) NxN matrices. The Haar measure ii{du^''^) 
for the unitary matrices -u*^'^ satisfies ^{du^^^) — Wg^ds^Xla^^jda^^^ds^^. Therefore the in- 
finitesimal variations of Y and 5* are given by 

/ i{dQV-dV Q)\ 

\i{dQV-dV Q) )^ 

( V 5u('^ V + Q 5u^''^ Q i (-V 5u^'^ Q+Q5u^^^ V)\ , . 

^[i{QSu^^^V-V5u'-''^ Q) QSu^'^ Q + V 5u^^^ V i' ^ ' 



dS = Y'^iSY - SY*)Y = Y'^{idM)Y . (3.12) 

This give us the real symmetric matrix dM in terms of the radial parameters {Aa}, the 
unitary matrices vM^ and their infinitesimal variations {d\a} and Su^''\ We just need to 
calculate a Jacobean which can be decomposed as the product of three determinants, 

n dMa,a+N = n ^r(^^^^ ' ^^-^^^^ 

a=l a=l ^V'^al-I- + 



N N 
,b+N 

a<b a<b 



n c?Ma,6+. c^M,,„+, = n 2 [J - ^ I c?a«dai? , (3.13b) 



AT AT I \ N / I \ I \ 



(3.13c) 



which eventually gives for the invariant measure of the symmetric unitary matrix 5" 

= a n - HM*'") (3.14) 

in terms of the measure iJt,{dX) — Y{a=i d^a of the matrix A and of the Haar measures //(du^'^) 
of the matrices w^'^ . 



IV. DYSON CIRCULAR ENSEMBLES AND QUASI-1£> DISORDERED 

SYSTEMS 

For the Dyson circular ensembles the number of ^'-matrices in a volume element dS of 
measure n{dS) around a given 5" is just proportional to n{dS) : 

P(dS) = ^ i,(dS) , (4.1) 
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is a normalization constant. Using the A-parametrization, one obtains that the matri- 
ces M*-'-* are independent from each other (except by symmetry relations) and distributed 
according the invariant Haar measure on the unitary group, while the N parameters 
are statistically independent from the w-matrices and have a joint probability distribution 
which can be expressed in the usual Coulomb gas analogy as a Gibbs function 

P({A4)=exp(-/37^({A4)) . (4.2) 

The symmetry parameter (3 plays the role of an inverse temperature and the effective hamil- 
tonian Ti. is characterized by a logarithmic pairwise interaction and a one-body potential: 

N N 

^({Aa}) = E /(Aa, A,) + E ^(Aa) , (4.3) 
a<b a=l 

f{Xa,Xb) = -log\Xa-Xb\ , (4.4) 



VW=(n + (^^ log(l + A). (4.5) 

This Coulomb gas analogy characterizes the orthogonal, unitary and symplectic ensem- 
bles which differ not only by the value of the "temperature" j3~^, but also by the presence of 
a small /3-dependent correction to the leading behavior of ^(A) in a large A^-expansion. It is 
remarkable that the pairwise interaction for the A-parameter is the same than in the global 



maximun entropy approach to the transfer matrix [15|, while V{X) differs in two important 
aspects: it is essentially proportional to log A instead of log^ A for large values of A, and 
the prefactor is just the number of modes instead of the classical conductance Nl/L^. 
These differences are not surprising since forward and backward scattering are essentially 
put on the same footing in the the circular ensembles, leading to a total transmission in- 
tensity of the order of the total reflection intensity, T R N/2 (up to weak-localization 



corrections) [p!7| , p!8| . For a disordered conductor or insulator, the refection R is much larger 
than the transmission T. Clearly, bulk diffusion characterized by an elastic mean-free-path 
/ cannnot be described by the circular ensembles, which are appropriate for systems where 
an injected carrier is subjected to a chaotic dynamics before finding [with equal chance] one 
of the two injection leads. However, as shown by Beenakker [^, in the large A^-limit, the 
A density-density correlation function (and therefore the variance of any linear statistics 
like the conductance) depends only on the pairwise interaction, and not on the particular 
form for V{X). Consequently, both the circular ensembles and the global maximum entropy 
approach to the transfer matrix yield identical UCF values 2/(16/3), slightly different 



from the perturbative microscopic result for quasi-ld disordered conductors 2/(15/?). 

In the polar representation of S", one can precisely see the difference between the circular 
ensembles and those appropriate for quasi-one dimensional disordered systems. For this, we 
just need to recall what we know from another statistical approach introduced for arbitrary 
N by Dorokhov from microscopic considerations and by Mello et al PE|PT[| from a 



maximum entropy assumption for the infinitesimal transfer matrix of the building block of 
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a quasi-l(i series. These works are based on an isotropy hypothesis: it is assumed that the 
M-matrices are distributed with the Haar measure on the unitary group and statistically 
independent from the radial part of M (see Eq. |2.8| ). This limits their conclusions to 
quasi-one dimension and yields a Fokker-Planck equation for P({Aa}) which implies the 
same UCF and weak-localization corrections for quasi-lrf conductors than those given by 
diagrammatic calculations. The evolution of P({Aa}) with the length Lx of the disordered 
part is given by an heat equation where the Laplacian becomes the radial part of the Laplace- 
Beltrami operator on a space of negative curvature. Using Sutherland's transformation, 
Beenakker and Rejaei have mapped this diffusion equation into a Schrodinger equation 
(with imaginary time) of a quantum set of point like particles free to move on a half line 
(the positive part of the real axis) within a certain potential. For arbitrary values of /3, 
these particles have a pairwise interaction, attractive for j3 = 1 and repulsive for /? = 4, 
making difficult to find the solution. Fortunately, this interaction vanishes for (3 = 2, and 
the solution of the diffusion equation is reduced to an exactly solvable quantum iV-body 
free fermion problem. This gives for the unitary case a pairwise interaction 



/(Aa - Afo) = - i In \Xa - h\ - ^ In 



arcsinh^ I a/A„ I — arcsinh^ ( \/A 



(4.6) 



which reduces to the usual logarithmic interaction assumed by the global approach for 
I Aa — Af,| <C 1, but which is halved if |Aa — Afe| S> 1 in the quasi-lrf diffusive or localized limit. 
This discrepancy is responsible for the slightly different UCF values characterizing ballistic 
quantum dots with chaotic dynamics and quasi-lci disordered conductors. 

In the localized regime, the global and local approaches give identical symmetry de- 
pendence of the localization lengths, though the (log) conductance fluctuations differs by a 
factor 2 in the quasi-ld localized limit [^. This later point again is consistent with the 
halving of the pairwise interaction /(Aa — Ab) for large eigenvalue separations given by the 
local approach. For metals and insulators far from a quasi-one dimensional shape, a more 



dramatic shrinkage of the validity of the universal RMT-correlations has been observed p4 



This means that transverse diffusion (or even more transverse localization |31]) yields a more 



signiflcant reduction of the RMT pairwise interaction than the one obtained in quasi-one 
dimension by Beenakker and Rejaei. To address this problem, we return to the study of the 
more familiar scattering phase shifts of S. 



V. SCATTERING PHASE SHIFTS IN QUASI II? METALS AND A SINGLE 

CIRCULAR ENSEMBLE. 

As discussed in the previous section, when applied to quasi-lc? conductors, the circu- 
lar ensembles give the right statistics for the u-matrices, and locally the correct interac- 
tion /(Aa,A{,), but certainly not the appropriate conflning potential ^(A). Using now the 
eigenvalue-eigenvector parametrization of S, we study the phase shift distribution in the 
case of weak disorder and quasi-lc? samples. We flrst introduce the basic notation, then 
write the expected universal correlation for Dyson ensembles, which we compare to our 
numerical results. 

Since the S'-matrix is unitary, its 2N eigenvalues are given by 2N phase shifts 6m- Fol- 
lowing Bliimel and Smilansky |]TU[ we write the phase shift density as 
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M 1 oo 

P{d) = E (^^^ - ^-)) = ^ E exp {-tne){TrS-) , (5.1) 

m=l 71=— oo 

where the angular brackets indicate average over the ensemble of disordered samples and 
M = 2N is the dimension of S. The two-point correlation function R2 is defined by 

M 

R2iei, 02) = E (^(^1 - (^^62 - dm')) . (5.2) 

When the phase shift distribution is uniform {p{6) = M/2ti) the two-point correlation 
function depends only on the difference rj = 62 — 61, 

M °° / 1 \ 



n=— 00 



The two-level cluster function is defined, for the reduced variable r = rjM/2'n', in the 
limit where the number of phases M goes to infinity, as: 

Y2{r) = lim ^/^(r) , (5.4) 

M^oo 



Using the expression (|5.3|) we have that 



>^'W^-^(l-2i:4'-(?=f^)) , (5.6) 



71=1 



where Fourier components are given by 

s: = ^{\TrS-\')-l. (5.7) 

The argument r of the cluster function goes from —00 to +00, and the Fourier transform 
of Y2, the two-level form factor (TLFF), is given by 

/oo 
dr Y2{r) exp {27rikr) . (5.8) 
-00 

Comparing the Fourier transform of Y2 and the Fourier coefficients of Y2 '^ in the large M 
limit, we can identify 

s*^ ^ - b{n/M) , M > 1 . (5.9) 

For matrices 5" belonging to Dyson ensembles, the distribution of phase shifts is given 
by the Coulomb gas analogy: 
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nR}) = ^exp(-/?m}) 



(5.10) 



where the effective hamiltonian is 



M 



^(R}) = -E log 

a<b 



(5.11) 



The cluster functions of the circular ensembles have universal forms which depends only 
on the system symmetries. For the unitary and orthogonal ensembles Y2 is an even function 
of its argument and has the form ^ 



sm vrr 



vrr 



(5.12a) 



,.nF/ X /sinvrrX^ ^^.^ ^ ^ /cos vrr sinvrrX s 
^2° W = - - ^e{r)) i^—^ - j^j . (5.12b) 

The corresponding form factors are also even functions of their argument, and have the 
universal forms: 

buE{k) = I J ^ if > 1 ' (5.13a) 



f 1 - 2fc + A; In (1 + 2fc) if A; < 1 

W/^) = I 1 _ ^ 1^ iffc>i • (5-13b) 

We check for a quasi-lci metal (Fig. 2) the agreement between the numerically generated 
Fourier components and the universal two- level form factors b{k) of Eq. (|5.13|) , assuming 
Eq. (|5.9|) . The S'-matrix of disordered strips described by a tight-binding Anderson model 
of 34 X 136 sites (details given in Appendix ^) are numerically evaluated. Averaging involves 
5000 different impurity configurations. The Fermi energy in units of the constant off-diagonal 
hopping term is E = —2.5, the number of propagating modes is = 14 and therefore the 
dimension of is M = 28. A relatively low wave- vector {k = 1.32 in Anderson units) is taken 
in order to avoid lattice effects, since we will be interested in the comparison between our 
simulations and analytical approaches (diagrammatic perturbation theory and semiclassical 
approximation) assuming a continuum limit. One finds a rather good agreement for the 
time-reversal symmetric case (run Rl, no magnetic field, COE-like) and for the non time- 
reversal symmetric case (run Fl, with magnetic field, CUE-like). The distribution of the 
phase shifts (inset) is quite uniform with and without magnetic field. The disorder in the 



^We follow the standard notation: Si(x) = dy; Ci{x) = C + Inx + /^f ^^^^ dy; e{x) = 

— 1/2, 0, 1/2, for X < 0, X = 0, X > respectively; C = 0.5772 ... is the Euler constant. 
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samples is very weak {W = 1 in units of the hopping term) giving an elastic mean-free-path 
/ = 0.7Ly, and an average conductance {g) = 4.1. 

We conclude that the phase shift density and correlations for a quasi-lrf conductor are 
well approximated by the corresponding COE-CUE density and correlations. Since it is 
clear in the A-parametrization that a quasi-lrf conductor cannot be seen as a member 
of a COE-CUE ensemble, we suspect that this numerical result merely indicates a good 
approximation, and that the main non COE-CUE behavior of S for a quasi-l(i conductor 
must occur in the distribution of its eigenvectors. As we will see in the following sections, 
the good agreement with the universal correlations gets poorer as we increase the disorder, 
enter into the quasi-lci localized regime, or go outside the quasi-lci geometry. 

VI. SCATTERING PHASE SHIFTS IN QUASI 1L> INSULATORS AND TWO 
UNCOUPLED CIRCULAR ENSEMBLES. 

The approximate COE-CUE behavior which we found in the previous section cannot 
remain in the presence of quasi-lc? localization for obvious reasons: -C 1 and a typical 
matrix element of the reflection matrices r and r' is much larger than those of t and t'. The 
matrix S can then be thought as two diagonal blocks, r and r', weakly coupled by t and t'. 
Using the polar decomposition, Eq. (pl9|), and the fact that the radial parameters {Xa} are 
exponentially large in the localized regime, one can write 

r' = ^(4)^(2) + 0(^-1) _ (Q lb) 

A strong quasi-lc? localization implies that S reduces to two uncoupled N x N unitary 
(symmetric if the presence of time reversal symmetry) matrices and isotropy means that each 
of them is invariant under orthogonal (unitary) transformation. The phase shifts associated 
to r and r' will then be described separately by two uncoupled COE-CUE ensembles. For 
a weaker quasi-lrf localization, we have a cross-over behavior between a set of 2N phase 
shifts with approximately COE-CUE correlations to two uncoupled sets of exactly COE- 
CUE phase shifts. We underhne that the observed COE-CUE phase shift distribution 
for the quasi-l(i conductor is less trivial than the COE-CUE character of the reflection 
matrices for strong quasi-lrf localization which only results from the isotropy assumption. A 
similar decoupling of the S'-matrix in two nearly independent blocks have also been discussed 
recently by Borgonovi and Guarneri |P3|| . 

If Y2{r) and b{k) are the two- level cluster function and the two- level form factor of two 
independent ensembles, the corresponding functions for the combined ensemble are given by 

m 

y^r) = \ y2(r/2) , (6.2) 
hs{k) = b{2k) . (6.3) 
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This expected crossover situation towards two uncoupled COE-CUE, for sample lengths 
of order of the localization length, is indeed observed in Fig. 3 where we show the Fourier 
components with and without magnetic field. For the very long, weakly disordered 
samples (RIO, diamonds, AR = 30, Vr=l, A^=14) the phase shift distribution (inset) is as 
homogeneous as for the quasi-ld conductor. The low harmonics of Y2 (small n values of 
s") behave as those of two uncoupled COE-CUE, reflecting the statistical independence 
of the short length trajectories contributing to reflection. Higher order harmonics behave 
more like a single COE-CUE, indicating that localization is not sufficient for decoupling the 
long trajectories. For a better understanding of this effect, one can use the semiclassical 
picture developed in Appendix y, where for relatively small values of n, is given in 
terms of periodic orbits of the closed sample which hit n times the vertical limiting hard 
walls: small values of n are related with short reflection trajectories which stay close to each 
of the samples edges and do not explore the other extreme of the sample (two decoupled 
COE-CUE behavior). On the contrary, for the large values of n, the contribution of the 
transmission trajectories can not be ignored, and removes the statistical independence of 
the two reflection matrices. 

Increasing of the system length being time consuming for the numerical simulations, 
we can more easily achieve localization by increasing the strenght of the disorder potential 
keeping the geometry fixed. In this case we are able to obtain a more complete decoupling, 
extending to higher harmonics the characteristic behavior of two independent COE-CUE 
(R4-F4, filled circles, AR = 4, W=4:, A^=14). Having stronger localization yielded by 
stronger disorder, we note an additional effect: the phase shift density is no longer uniform. 
We will discuss this departure from isotropy in detail in the next section. For the purpose 
of the present discussion we only indicate that we numerically unfold the phase shift spectra 
to a rescaled spectra of uniform density. Then, one can see in Fig. 3 a rather complete 
statistical decoupling of the two sample edges. Let us note also that the small n behavior 
of s^' is now a little above what we expect from two uncoupled COE-CUE, a point which 
will be considered in Sec. |VII1| . 

In order to study more precisely how the transition from the metallic case (COE or 
CUE-like cases) to the localized regime takes place, we calculate the number n{r) of levels 
contained in an interval of length r for the unfolded spectum of the phase shifts and its 
variance (number variance) 

S2(r) = ((n(r) -r)^) , (6.4) 

which can be obtained directly from the numerical data and compared with universal forms 
of the standard ensembles. Using Eq. ( p.l2| ), one gets for the unitary and orthogonal cases 



Sf/F(^) = -t: (In (2Tcr) + C + 1 — cos (27rr) — Ci(27rr)) + r(l- -Si(27rr)^ , (6.5a) 



E|,W = 2E?„W+f^) ,6.5b) 

\ TT / TT 



with the large-r behavior 
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(6.6a) 



Soe(^) = ^ (in (2vrr) + C + 1 - n^s) + 0{r-') 



(6.6b) 



For the superposition of two independent ensembles, one gets |^| 



S^(r) = 2 S2(r/2) . 



(6.7) 



Since our original phase shift spectrum is bounded between and 27r, S^(r) folds back to 
for r = M as the number of phases in < < 27r (or < O < M) is always M. Hence our 
comparisons between our numerical data and Eqs. (|6.5| ) - ( |6.71 ) are meaningful for r <^ M 
and we focus our attention to intervals < r < M/3. 

In Fig. 4 we show how the number variance S^(r) changes from the metallic to the 
localized regime when we increase the disorder or the sample length. If the disorder is 
increased, we obtain a continuous transition from the COE case towards less rigid spectrums 
and the phase shift density develops a more and more non uniform structure. The strength 
of the local disorder not only decouples the right and left reflections, but also introduces in 
the transverse direction a dimensionality dependent dynamics, which breaks isotropy. The 
number variance can then exceed that of two uncoupled COE. 

When the disorder is small, and localization is achieved by increasing the system length, 
the transverse dynamics remains essentially zero dimensional, S is isotropic and the number 
variance indicates a cross-over from approximately a single COE for the 2N phase shifts 
of 5* (quasi-lci conductor) to two decoupled COE-like sets of N phase shifts associated to 
right and left reflections. 

When a magnetic field is applied, we obtain similar results, with a slightly improved 
agreement with the CUE-like character: the applied magnetic field suppressing coherent 
interferences between time reversed trajectories shifts the sample towards the metallic regime 
and doubles the localization length 

VII. AVERAGE SCATTERING MATRIX AND ISOTROPY HYPOTHESIS 

Quasi-lc? distributions are partly based on the isotropy hypothesis: i.e. the unitary 
matrices m*^'-* are uncorrelated with the radial part of S and distributed with the invariant 
measure on the unitary group. This implies that the ensemble average of S must be zero. 
This property does not hold for high disorder and non quasi-lci geometry. For instance, a 
non uniformity of the phase shift distribution was noticed in the previous section for a large 
disorder and occurs for weaker disorder in shorter samples (squares and thin slabs). This 
source of discrepancy with the circular ensemble is studied in this section. 

In the circular ensembles, the phase shift density is uniform and equal to po = M/27r. 
This means that (TrS"") = for all n ^ 0. We will see that this is not the case for disordered 
systems. The first harmonic of the Fourier expansion of the phase shift density p{6) is: 



1 



Re[exp{-i9){TrS)] , 



(7.1) 
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and (TrS) = is a necessary condition to obtain uniform phase shift distribution. Since 



N 



{TrS) = Y.i{raa) + {r'J), (7.2) 



a=l 



we can calculate pi(^) by evaluating the mean values of the diagonal reflection elements in 
perturbation theory. 

The transmission (reflection) amplitude from a mode a on the left to a mode b on the 
right (left) for electrons at the Fermi energy Ej = h'^k'^/2m is given by 



ha = -iTiiyaVhf^ / dy' / dy (pliy') (f)a(y) Gk{L^,y']0,y) (7.3a) 



Ha = Sa,b -ih{vaVhY^'^ / dy' / dy (pHy) (j)a{y) Gk{0,y';0,y) (7.3b) 



where Va (f f,) and (pa {4>b) are the longitudinal velocity and transverse wavefunction for the 
incoming (outgoing) mode a (b). For hard- wall boundary conditions, the transverse wave 
functions have the sinusoidal form presented in Sec. |I|, Vn = hkn/m, k"^ = k"^ — (nn/Ly)^, 
n = a,b. We note by m the effective mass of the electrons. For the transmission (reflection) 
amplitudes Gk{r'; r) is the retarded Green function evaluated at the Fermi energy between 
points r = [x, y) on the left lead and r' = {x' , y') on the right (left) lead. Similar expressions 
hold for the transmission (reflection) amplitudes for modes comming from the right by using 
Gk{0, y'] L^, y) in Eq. ( |7.3a| ) and Gk{L^, y'; L^, y) in Eq. ( [7.3 b| ) instead of Gk{L^, y'; 0, y) and 
Gk{0, y'] 0, y) (and placing the y abscissa at a; = L^., and y' at x = {L^)). 

For a given impurity configuration, the unaveraged retarded Green function Gk for elec- 
trons at the Fermi level in the absence of a magnetic field satisfies 

(1^ V? +^ + Vir) + .7) Gk{r'; r) = 5(r' - r) , (7.4) 

with 7 — > 0"*". We will assume that the impurity potential V{r) is given by Nt uncorrelated 
(5-function scatterers (of strength u) randomly distributed in the disordered strip. 

V{r) = f:u6{r-R^) {I'^ff'^^i; « = !,.. .,iV,. (7.5) 

— 1 \ ^'iV y 

The standard techinque used in disordered systems is to solve Eq. ( [7.4| ) in perturbation 
theory and take the ensemble average at each order of the perturbation expansion |25|,0. As 
indicated diagramatically in Fig. 5, we merely expand to second order in the perturbation 
expansion. 

The unperturbed Green function G^j^\r'; r) for an infinite strip can be expanded in the 



base of the transverse wavefunctions and the n^^ coefficient is a one- dimensional Green 



function with an effective wave vector /c„ = y/c^ — [mi/LyY {kn is real when n corresponds 
to a propagating mode and pure imaginary otherwise). 
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.^2 H TT - a;|)0„(y')0n(y) 



(7.6) 



The inclusion of G^^ in (|7.3b|) just cancels the 5^,6 factor since in the absence of disorder 
the modes propagate without any reflection. The first and second order corrections are 
respectively given by 



(7.7a) 



where Ui = Ni/{LxLy) is the impurity density. For the second order correction we are 
only giving the leading term (in the impurity parameter urii, and in the inverse mode 
number 1/A^), and we have cut the sum over the internal momentum (as usually done for 
5-function potentials). The highly oscillating phases e^*'^"^=" would be suppressed by inelastic 
scattering. For translationally invariant (after averaging) systems the perturbation theory 
is usually done in momentum representation; the first order term gives rise to a real self- 
energy that merely renormalizes the Fermi energy, while the second order term gives rise 
to an imaginary self energy which is responsible for an exponential damping of the average 
real space Green function. For the reflection amplitudes from finite disordered regions, we 
take into account perturbation up to second order term only. This calculation neglects 
the multiple scattering processes characteristic of the diffusive regime. However, we are only 
interested in the average value of the reflection amplitude, where single scattering dominates 
multiple scattering, as we will see in our semiclassical approach, and Eqs. ( [TTI ) give the 
main features seen in the numerical simulations. They are proportional to the strength of 
the disorder and show a non-zero average only for the diagonal elements of the S'-matrix. 
The denominator ka indicates that the absence of self-averaging is more pronounced for the 
higher modes. The breakdown of perturbation theory for small ka is understandable since 
it is close to a threshold of complete reflection which cannot be obtained perturbatively. 
Clearly, those threshold effects, when new conduction channels appear (A^ —>■ + 1)), 
drastically limit the validity of the isotropy assumption. 

Up to second order perturbation the average diagonal reflection amplitudes are 

= l«l ^ V2(l-cos(2fc.L.)) , (7.8) 

where a = (^))- For the lowest modes Eq. (fTSp gives a correction vanishing 

as l/k"^ (or {Ly/NY), but the correction remains important for the highest modes. On the 
other hand, in our numerical simulations in a lattice model, will always be finite and 
not very large. In Fig. 6. a we show the values of K^aa)! obtained from our numerical 
simulations, as a function of the mode number a. The average reflection amplitudes of 
the sample described in the Sec. (Rl, squares, AR = 4, W=l, N=14) are close to the 
functional form 1/ka (solid thick line) and to those of a longer sample (R8, circles, AR = 10, 
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W=l, N=14). The average value of the nondiagonal reflection amplitudes (not shown) are 
zero within the statistical error. Since this perturbation calculation is performed in the 
continuum and yields a rapidly oscillating phase we do not expect to get full agreement with 
the numerical simulations on a lattice. We are mainly interested in the mean behavior of the 
diagonal reflection elements as a function of the channel number a, the number of modes 
and the strenght of the disorder W. Increasing the disorder (R2, filled diamonds, AR = 4, 
W = 1.5, = 14) enhances the correction. Keeping the number of modes = 14 (Fig. 
6.b) while going to a square geometry (R12, squares, AR = 1, = 1, = 14) decreases the 
Fermi momentum and gives a larger correction for the average reflection amplitude than for 
the sample of Sec. |V[ The correction for a square geometry is of the same order of that of 
a long sample (Rll, circles, AR = 4, W = 1, N = 14) with the same number of modes and 
cross section (same Fermi momentum). Increasing the disorder (R13, diamonds, AR = 1, 
W = 1.5, A^ = 14) enhances the correction, while increasing the number of modes (R18, full 
triangles, AR = 1, W = 1.5, A^ = 28) decreases the correction. 

The mean values of the scattering matrix elements have been considered by lida, Wei- 
denmiiller and Zuk who studied the interplay between universal conductance fluctuations 
and the statistical properties of the hamiltonian spectrum of a disordered strip. They define 
the sticking probability Ra = 1— |(?"aa)P; measuring the weight of the fast processes (where 
the particle is reflected after a few scattering events in the disordered region) versus the long 
trajectories. Sticking probabilities smaller than 1 give corrections to the universal quasi-lD 
value of the variance of the conductance fluctuations, which can be calculated with the aid 
ofEq. (|73). 

This non-zero average of the diagonal refection elements is in disagreement with the stan- 
dard maximum entropy approaches where one assumes that the matrices u^^^ are uniformely 
distributed in the unitary group, giving zero average values of the S'-matrix. Recent work by 
Mello and Tomsovic has relaxed the isotropy assumption, making possible a non-zero 
mean for the transmission amplitude (diagonal in mode number and exponential in l/Lx). 
Our results show that a non-zero average for the diagonal reflection amplitudes is needed in 
order to describe systems outside the weak disorder and quasi-one- dimension. 

In the large A^-limit we approximate the second order perturbation calculation of (TrS) 
by converting the sum over modes into an integral, the phases exp {2ikaLx) give rise to 
higher order terms in 1/A^, and from Eqs. ( [7.7| ) we obtain (away from the thresholds) 

(Tr5(^'^)) ^ a ^ In f P^^^] - « ^ ^ In (2Ar) . (7.9) 

It is important to notice that Eq. (|7.9|) does not depend on the length of the sample 
(or the aspect ratio) but only depends on the properties at the entrance of the sample 
(transverse cross-section Ly and Fermi momentum k). Its logarithmic dependence on A^ 
indicates that there is homogeneity in the phase distribution only for very large A^ and very 
small transverse cross section. 

From Eqs. ( [7. 1|) and (|7.9| ) we can see that the phase shift distribution will not be uniform, 
unless we have a narrow sample with weak disorder and a large number of propagating modes. 
A non uniform density is often the case for hamiltonian spectra (or for the radial parameters 
of M), and one studies the correlations of the unfolded spectra JTT|J5D[] . Instead of studying 
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the bare phase shifts 9m, we have to consider the rescaled variable 9m = M{9m), where 
AfiO) is the number of levels below 6, i.e. AfiO) = p{6')d6'. Since our interval of phase 
shifts is bounded between and 27r, for the numerical study, we repeat twice each phase 
shift sequence and calculate the two-point correlation function in the second interval (M to 
2M for the unfolded variables Qm) ■ The Fourier components can then be compared with 
the universal two- level form factors of Eq. ( ^.13|) , as we do in the following sections. 



VIII. UNFOLDED PHASE SHIFT SPECTRA OUTSIDE 
QUASI-ONE-DIMENSION 

We now turn our attention to samples far from a quasi-lci shape. For disordered squares 
(Fig. 7. a) the phase shift densities are non uniform in the absence of magnetic field (inset). 



The degree of nonuniformity is consistent with Eqs. (|7. 1| ) and (|7.9|) , i.e. it becomes more 
important for higher disorder and smaller number of propagating channels. The Fourier 
components (after unfolding) are in relatively good agreement with the COE-CUE two- 
level form factors, but the correspondence becomes poorer when we increase the disorder. 
For the large-n Fourier components, we are probing the correlations for small separations, 
where we can take the phase shift density as constant. This approximate translational 
invariance makes that only the term with n = —n' survives in the expansion of the delta 
functions of Eq. ( |5.2|) . The Fourier components of the two-point correlation function are 
still given by (iTrS*"!^), through Eq. ( |5.7| ). For squares with magnetic field and low disorder 
the phase shift density is relatively uniform (lower histogram of the inset) and we use Eq. 
( p.7[ ) instead of unfolding the spectrum. 

In Fig. 7.b, a disordered slab is considered (upper inset, the incoming electron flux is along 
the direction of the shortest dimension). The phase shift distributions (inset) are strongly 
non uniform {Ly in Eq. ( [7.9|) is very large), but the unfolded spectra, with (triangles) and 
without (filled squares) magnetic field, are relatively well described by Eq. ( ^.13|) , except 
for small n values where a careful look indicates values above the COE-CUE behavior. This 
is another important source of departure from the universal behavior related to transverse 
diffusion, in addition to the cross-over mentioned in section |VI| coming from longitudinal 
localization. 

The approximate agreement of the numerical data with the two-point form factors of 
the circular ensembles is somewhat surprising in these geometries. Given the approximate 
agreement of the unfolded spectra with the standard ensembles at the level of the two- 
point correlation function, we might ask at this stage whether the unfolded phase shift 
spectrum of metallic conductors is well described by the circular ensembles, independently 
of the shape and strenght of the disorder. However, a check at the level of the two-point 
correlation function is not very accurate, as we have learned from statistical studies of chaotic 
hamiltonians [111] and transmission matrices p3. A better test is provided by integrals 



involving the two-point correlation function. As in Sec. |V|, we now consider the number 
statistics n{r) and the number variance (Eq. ( |6.4| )) 

In order to analyze systematic departures from the random matrix correlations we plot 
in Fig. 8, for various geometries, degrees of disorder and number of propagating modes, 
the difference of the number variances between the numerical data and the COE-CUE 
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values ( |6.5|) . We consider a^{r) = S^(r) — T,Q^{r) for samples without magnetic field and 
(T^(r) = S^(r) — E^^(r) for samples with nonzero magnetic field. Like in the previous section, 
we only show the interval < r < M/3 where the comparission is meaningful. One can see 
very clearly now that, even after unfolding, the accuracy of the random matrix description 
strongly depends on the shape and degree of disorder of the samples. The difference (j^(r) 
in number variances grows approximate linearly with r. The magnitude of the slope of cr^(r) 
measures the validity of the random matrix description. Squares with low disorder (after 
unfolding) are well represented by the standard ensembles, the agreement becomes poorer 
when increasing the disorder, and improves when augmenting the number of modes A^. 
However, slab-shaped samples show also large deviations respect to the COE-CUE values. 
The conditions for having a good COE-CUE distribution coincide with the conditions for 
diminishing | (TrS) \ and obtaining a more uniform distribution (Eqs. ( [7.1| ) and ( |7.9| )). Notice 
however that Fig. 8 is done over unfolded spectra, where the nonhomogeneous phase shift 
distribution is in principle already accounted for. Our task in the remaining section is to to 
quantitatively study the deviations to the random matrix behaviors in disordered conductors, 
giving the slope of the curves a^{r) as a function of the aspect ratio, degree of disorder and 
number of propagating modes of the disordered sample. 



IX. DEVIATIONS FROM THE UNIVERSAL RMT BEHAVIOR 

The statistical properties of the spectra of small disordered systems have been consid- 
erably more studied than those of the scattering phase shifts. It is then appropriate to 
establish a connection between the two in order to understand our numerical results. The 
relationship between phase shifts and energy levels has been investigated in the semiclassical 
limit by Bogomolny |^ and by Doron and Smilansky . Their approach is not completely 



applicable to our case but provides some guiding concepts. We reproduce in Appendix ^ 
their main results and discuss the points where the correspondence does not hold. 

Under the assumptions discussed in Appendix ^ of neglecting the evanescent modes and 
the nonuniformity of the phase shift distribution, the energy level density is given by 

d{E) = ^{t^{E)) , (9.1) 

where Tw{E) is the Wigner time 

r.,E,.±Tr[sW-^), (0.2) 

whose physical interpretation is discussed in Appendices P and Assuming that the 
density of states of our disordered rectangular samples is the same as without disordered 
d{E) = 1/A = mLxLy/27rh'^, we can check the validity of Eq. (|9.1|) ) against our numerical 
simulations. As indicated in Table I of Appendix ^ there is agreement within 5 to 25 % in 
most of the cases, which is reasonable given the various approximations involved. 

The quantity most often calculated in metallic spectra is the density-density correlation 
function 
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K2{E,E + e) = Y^{5{E-En)5{E + e-En')) - d{E) d{E + e) (9.3) 

n,n' 

Scaling the energy separation with the mean level spacing (e = e(\{E)) and the phase 
shift separation with the mean phase shift distance (r = rjM /2tx) we have, from Eq. ( |B8D , 
in the large M limit that 

K^{e/d{E))d-\E) = 5{r) - K,(r) , (9.4) 

As stated before, we will put aside the fact that the two above mentioned assumptions 
are not quite true in our systems and we will pursue the consequences of Eq. ( |9.4| ) . 

The density-density correlation function for disordered systems has been obtained in 
perturbation theory by Altshuler and Shklovskii [Q 

K2{e) = -Rey^{e + i?iDq^ + (9.5) 

for energies e large compared to the level spacing A, and small compared with the energy 
scale h/Te, associated with the elastic scattering time Te- The factor s accounts for the spin 
degeneracy of each level (s = 1 since we work with spinless electrons), and 7 is a small 
energy cutoff (to account for inelastic scattering). For simplicity we will be work the case 
of zero magnetic field. The sum is over the diffusion modes in the sample, assumed to be 
a dimensional parallelepiped with sides L^, that is, = 'rc'^J2'l=i{n^/ L^y. The diffusion 
coefficient is D = Vfl/d. The Thouless energy E^.^^ = hD/L'j^ is inversely proportional to the 
time that takes an electron to diffuse across the sample in the /x-direction. For samples with 
all equal (hypercubes) we just have one Thouless energy. That will be the case of our 
square samples, while for quasi-one-dimensional samples we have in principle two Thouless 
energies, but we will reserve this name for the smaller one, that is, the one associated with 
the length L^. 

The mean square fluctuation in the number of levels (number variance) is given in terms 
of the density- density correlation function, 

fE+e/2 fE+e/2 

{[5N{e)f) = dEj dE2 K2{Ei,E2) , (9.6) 

JE-e/2 JE-e/2 

which from ( p.5| ) can be written as 

m{e)f) = ^ E In 

For energies e <^ E^^^ the sum ( |9.7[ ) is dominated by the term with all null and [^] 

([5Ar(.)f) = ^ln(£^ + l) . (9.8) 

The perturbation theory of Ref. is valid for energy separations e larger than the 
inelastic scattering 7 or the level spacing A. In our case there is no inelastic scattering 



(7 + fiDq^f 



+ 1 



(9.7) 



22 



and we substitute 7 by A obtaining the standard random matrix theory result {[SN{e)Y) = 
2 / ir'^ In {e / A) for A <^ e <^ Ec- For energies e ^ the summation over {n^} can be 
replaced by an integral over (ig^ and 

([5iV(^)f) = c.(-^) , (9.9) 

where Cd is a dimensionality-dependent numerical coefficient The asymptotic results 
( |9.8|) and (|9.9|) have recently been rederived by Argaman, Imry and Smilansky using 



a more intuitive semiclassical method for electrons in the diffusive regime. The small- 
energy universal regime is obtained for times long enough to allow a diffusing electron 
to ergodically explore all the sample, while for energy intervals larger than the Thouless 
energy (short times) the correlations depend on the diffusing, unbounded dynamics of the 
electron, which is dimensionality and disorder dependent. The universal character of the 
short range eigenvalue correlation and the long range non-universal (dimensionality and 
disorder dependent) part have been also obtained in numerical simulations of by Dupuis 
and Montambaux ||^ who studied the crossover between the two regimes in an Anderson 
model. 



Given Eq. (O), and the above results for the number variance predicted by perturbation 
theory, we expect the correlation functions of the phase shifts to have a COE-CUE behavior 
for phase shift separations smaller than 

Vc="-^. (9.10) 

Our numerical data summarized in Table I of Appendix A indicates that for the studied 
samples the critical angle r]c is smaller than the mean phase shift spacing {tc = r]c/ (27r/M) = 
{g) /2tt < 1). Therefore they are in the transition regime between the two asymptotic limits 
(PD and (lop. 

In Sees. and |VIII| we have studied the number variance of the phase shifts cr^(r) = 
E^(r) — IlQ^{r). In particular, from Fig. 8 we can see that cr^(r) grows almost linearly 
from r ^ r^. To check our prediction, it is interesting to look at the value of the slope. 
Whithin the perturbation theory of Ref. Q the difference between the number variance for 
the real energy levels and the R.M.T variance is obtained by excluding the term = riy = 
from the sum ( |9.7] ). For a two-dimensional sample the slope of the difference evaluated at 
z/c = Ec/ A is given by 

For a quasi-one- dimensional sample the sum is over ^ and the slope is approximately 
0.0045/ec. Given the relationship between the statistics of phase shifts and energy levels, we 
expect that the slope (f)c = ^o'^('^)|r=rc ^-Iso scales linearly with l/vc- In Fig. 9 we show the 
value 0c versus the inverse Thouless energy l/vc- The linear relationship predicted by Eq. 
( p.ll| ) turns out to be approximately valid. The slope obtained for the two-dimensional case 
(squares geometries in Fig. 9) agrees within 50% with the coefficient of ( |9.11|) and is a factor 
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of 4 larger than the slope obtained for quasi-one dimensional geometries (roughly the same 
ratio than for the eigenenergies) . Obviously, in those samples, the quasi-one-dimensional 
limit is not achieved and the approach to this limit depends on the aspect ratio of the 
sample. 



X. CONCLUSIONS 

In this work we have considered scattering in disordered systems. We calculated the in- 
variant measure of S in the polar decomposition and compared the COE-CUE distributions 
with those of the global and local approaches to the transfer matrix of disordered systems. 
We then turned to the study of the scattering phase shifts for different geometries and de- 
gree of disorder. For quasi-one dimensional samples in the metallic regime the density and 
the two-point correlation functions are close to those of the circular ensembles. Increasing 
the lenght with fixed disorder or the disorder with fixed geometry, we break 5* into two 
uncorrelated reflection submatrices. Deep in the quasi-lc? localized regime the decoupling 
is almost complete and the two-point correlation function of the phase shifts is that of a 
superposition of two independent circular ensembles (COE or CUE). 

Outside quasi-one dimension, the phase shift density differs from the uniform distribu- 
tion characteristic of Dyson ensembles. This anisotropy comes from short-time processes 
which yield a non vanishing ensemble average of S, and is strongly enhanced in the vicinity 
of a energy threshold where a new conduction channel appears. Using diagrammatic per- 
turbation theory, we calculated the average value of the reflection amplitudes, finding good 
agreement with the numerical simulations and establishing the conditions for having the 
uniform phase shift distribution of the circular ensembles. One needs both energy far from 
a threshold, weak-disorder, small transverse length and large number of propagating modes. 
When this is not the case, anisotropy is important and we have to unfold the spectrum in 
order to study its correlations. 

For large transverse lenghts the electron diffusion in the perpendicular direction gives rise 
to small-fc corrections of the two-level from factor b{k) of the [unfolded] phase shift spectrum. 
This non universal behavior is related to the large-energy (non-ergodic) behavior found by 
Altshuler and Shklovskii for the energy-level statistics of small metallic particles. This lead 
us to relate under certain assumptions the scattering phase shift and the energy-level form 
factors. We numerically checked this relationship for the number variance (fluctuation of the 
number of states within a given interval). We verified that the non universal discrepancies 
become increasingly important as the Thouless energy decreases, in a geometry dependent 
fashion. 

One of the most important problems remaining in the random matrix theory of electron 
transport in disorder systems is the extension of the standard approach outside the quasi- 
one-dimensional case. We have shown in this work that anisotropy is an essential ingredient 
in higher dimensions. The precise form of the decoupling of the S matrix into independent 
reflections blocks in the one- dimensional localized regime is another interesting problem 
left for futures studies. Given the recent progress in the theory of parametric correlations 
1^-^ , the energy dependence of the scattering phase shifts desserves further studies (e. 



the correlations of the characteristic times t^^E) (Eq. (p^), slopes in Fig. 10) which we are 
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presently developing. 
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APPENDIX A: SUMMARY OF THE NUMERICAL SIMULATIONS 

In this appendix we summarize the results of the numerical simulations and we briefly 
indicate the way they were obtained. The simulations were performed in a tight-binding 
model with a number of sites x Ly and a random on-site disorder of amplitude W by using 



a recursive Geen function method P7| , ^ . In a magnetic field we use the Peierls substitution 
to relate the hopping matrix element to the vector potential. The magnetic field is taken to 
be linearly increasing in the leads from zero to its full value in the disordered region over a 
distance of the order of the transverse dimension Ly. For the small fields that we work with 
the statistical results are independent on the way the field is introduced. 

Runs 1-11 in Table I are in the quasi-one- dimensional limit, while runs 12-18 are for 
square samples and runs 19 for thin slabs. The time limitation for the numerical simulations 
is the transverse size and that limits the number of samples {NS) considered for large Ly. 
The number of modes was kept low (14-20 in most of the cases) in order to have small 
Fermi energy {Ef) and wavevector {k) avoiding lattice effects. The Wigner time r^- is 
calculated from Eq. ( |9.2[ ) and provides, through (|9.1|), the value of the mean level spacing. 
Using the free-space two-dimensional density of states to give the Weyl term of our disorder 
region, the level spacing is A = mLxLy/27ih'^, and the column showing [ry^^/K) (N/tt) 
checks the approximate validity of Eq. ( p.l| ) under the hypothesis that the level spacing 
does not change by the effect of confinement or the disorder. We obtain in most of the cases 
an agreement within 5-25 % is quite reasonable given the various approximations involved. 

The elastic-mean-free path / calculated in Born approximation using samples of 
transverse dimension Ly and a few lattice sites on the longitudinal dimension. In the metallic 
regime it agrees very well with the mean free path obtained from the Drude formula for the 
dimensionless conductance (g) = NItt/2Lx. The Thouless energy in Anderson units is given 
hy Ec = 1/ Ll. sinfc. The reduced critical angle Tc is obtained by scaling the critical angle 
r]c = TwEjh (Eq. ( CT) ) with the mean level spacing 2tt/M. The factor (pc measuring the 
discrepancy with the standard ensembles is obtained from Fig. 8 by taking the slope of (T^(r) 
between r = 0.7 and r = 7. This is not completely equivalent to the slope at since Tc and 
the range of linear behavior of cr^ varies from sample to sample, but has the advantage of 
providing a consistent definition for all the cases. 

Conductances for individual samples are obtained from Landauer's formula, Eq. ( p.. 21 ), 
and then averaged over impurity realizations yielding the average {g) and the variance {Sg"^). 
According to the Thouless formula {g) = 27r£'c/A in the metallic regime, and by using 
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the approximate A of the perfect square we obtain an agreement of 25% with Landauer's 
formula. It is important to notice that is obtained from / and r^y, directly extracted 
from the simulations (without any assumption over the form of the conductance). In the 
localized regime (runs 4,5 and 10) the average conductance and its variance are dominated 
by special configurations with anomalously large values. Therefore averaging over log is the 
meaningfull thing to do. In Table I we give as {g) for these runs the typical conductance 
exp (log (7) and {5{\oggY) in the column of the variances. 

In the metallic regime the difference between the average conductance without and with 
magnetic field is the weak locatization correction. Mello and Stone |^D[ have calculated this 
correction in the quasi-ld limit by diagrammatic and random matrix approaches and found a 
value of 1/3, which is somehow larger than the one obtained from Table I. It is rather difficult 
to reproduce the theoretical value in numerical simulations on disordered strips since getting 
to the regime of validity of the theory (diffusive, quasi-lD, metallic and magnetic field large 
enough to kill the time reversal symmetry but small enough for not to alter the classical 



electron paths) requires extremely large sample sizes ^T|- The values of the conductance 
fiuctuations {^g"^) for zero magnetic field in the metallic regime are in good agreement with 



the diagrammatic values and previous numerical simulations p],^. The halfing of the 



variance by the effect of the magnetic field is not quite obtained for 0/0o = 2, but for ratios 
of the order of 5-10 the factor of 2 it is always found. This is somehow surprising since for 
0/00 = 2 we have already a complete transition towards the unitary ensemble (see Fig. 2). 

The mean values of | (TrS) \ are in qualitative agreement with the perturbation result of 
Eq. (|7.9|) , increasing with disorder and decreasing with the number of modes N . 



APPENDIX B: RELATION BETWEEN THE SCATTERING PHASE SHIFT AND 

ENERGY LEVEL STATISTICS 

In this appendix we present, following Bogomolny and Doron and Smilansky 



a relation between the statistical correlations of the phase shifts and the energy levels in a 
form appropriate to the scattering which we consider. 

Close to the entrance cross sections (x < and x > L^) the scattering wave functions 
are given by the asymptotic expressions ( p.l|) plus the contribution of the evanescent modes: 



A'' 1 00 -| 



y) = T.TTH (^ne^'"^ + i^ne-^'^"^) Uv) + E urn7^^"el'^"l>„(t/) , (Bla) 

n=l kn n=N+l l'*^"! 



Af 1 00 1 

n=l kn n=N+l ^^n] 

(Bib) 



In Ref. [^, the contribution of the evanescent modes for the evaluation of \l//(x = 0, y) 
and n{x = Lx,y) is neglected in the semiclassical limit. This sharp cut off at + 1 
in the sum over n is clearly not exact in the case of our numerical simulations, where 

is typically 14 and higher transverse harmonics are needed in order to represent the 
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scattering wavefunction at the entrance. Neglecting the evanescent channels, the scattering 
wave function vanishes at the boundaries (x = and x = L^) ii B = —A and C = —D, 
implying that the scattering wave function coincides inside the sample with an eigenfunction 
of the hamiltonian (assuming hard-wall boundary conditions). Therefore, an energy E is an 



eigenvalue of the closed system whenever S{E) has an eigenvalue — 1 pT 



det{I + S{E)) = . (B2) 

As illustrated in Fig. 10, the energy levels can be obtained by the intersection of the curves 
9m{E) with horizontal lines at 6' = (2n+l)7r. The accuracy of this quantization condition has 
been recently tested numerically in the Sinai billiard. Away from the mode thresholds, 
the errors in the determination of the energy levels are of the order of only a few percent, 
and they are greatly reduced when a few closed channels are kept. On the other hand, we 
are not interested in the precise energy levels but in the relationships between their statistics 
and that of the phase shifts. Towards this end, we have calculated the nearest neighbour 
distribution of the energy levels determined from ( P2|) and found, in the metallic regime, 
good agreement with the corresponding Wigner distributions. From the above quantization 
condition (p2|) Q the level density can be expressed as 

1 M oo 

m = -E^)) = —-Y: E i-^TME) exp {in9^{E))) , (B3) 

n m=l n=-oo 

where the time associated with the m}^ phase shift is 

rJE) = n . (B4) 

Ignoring correlations between 6m and and in the hypothesis that the phase shift 
distribution is uniform {{TrS^) = for all n ^ 0) we have Eq. of the text 

d{E) = ^ ME)) , (B5) 
where the Wigner time Tyy{E) is given by ( pl2| ) or 

1 M 

^^-(E) = TT E ^m{E) . (B6) 
m=l 



If we use the convention 5* = and we ignore the evanescent modes, the condition 

\r t) 

to have an eigenenergy of the closed system with periodic boundary conditions is that S has the 

eigenvalue 1, that is, det(I — S{E)) = 0. It is easy to see that S = SU, with U = 

Therefore, for Dyson Circular Ensembles where the probability distribution is given by the invariant 
measure in the appropriate space, we do not expect a difference between the statistical properties 
of S and S. 
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The Wigner time is the mean slope of the curves 9m{E) in Fig. 10 and has the physical 
interpretation of the mean time spent by the particles in the scattering region, as can be easily 
seen from our semiclassical analysis of the transmission amplitude (Appendix |C|). Notice 
that the Wigner time resulting from our convention for the S matrix is not the usual delay 
time but and absolute (always positive) time. In particular, for a perfect (non-disordered) 
sample it is not zero, but given by = 2h/MJ2n=id{knLx)/dE = jfY^n=i^x/vn, that is, 
the average over the propagating modes of the ballistic traversal times. 

From its definition ( |9.3|) and the relationship ( p3|) we can write the density-density 
correlation function of the energy levels as 

1 M oo 

' m,m'=l n,n'=— oo 

{t^{E)t,^,{E) exp {i{eUE)n - 9^,{E)n' - r„,{E)n' e / Ti)) - <\{E)d{E + e) , (B7) 

where we have used 6m' {E + e) ^ 9mi{E) + Tm'{E)e/h for e <^ E and we assumed that 
the Tm{E) are smooth functions (in the scale of e) which can he taken outside the ensemble 
average. If the phase shifts are uniformely distributed only the terms with n = n' survive 
the ensemble average and [^0| , ^ 



mm + iHE) (^j^fi. (|^d(j5)) - 1) . (B8) 



Scaling the energy difference with the mean level spacing (e = ed{E)), in the large 
M limit we obtain ( |9.4|) of the text implying that the phase shifts and the energy levels 
have the same two-point correlation functions when expressed in terms of their respective 
reduced variables (r and e) . Notice that the relationship between the phase shift correlations 
(given by R2{r) or Y2{r)) and the energy level correlations (given by i^2(e)) is obtained in 
two steps: first the quantization criterion ( |B2|) (which neglects the evanescent modes) and 
second assumptions on the homogeneity of the phase shift distribution and the smooth 
character of Tm{E). 



APPENDIX C: SEMICLASSICAL APPROACH 

In this appendix we present a semiclassical treatment of the scattering in a wave-guide 
geometry with disorder, for understanding of the non-zero average of the diagonal reflection 
amplitudes and of the various time scales in the problem. The semiclassical approximation to 
the transmission and reflection amplitudes can be obtained by replacing the Green functions 
Gk of Eqs. (|7.3| ) by their semiclassical path-integral expression ||5^ 

^ ' s{y,y') 

The sum is over classical trajectories s between the arguments of the Green function (we 
will be omiting the x dependence since it will be assumed x = or L^, depending whether 
we are considering reflection or transmission), Ss is the action integral along the path s at 
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energy Ef = h^k'^ /2m = mv'^/2, the stability prefactor is Dg = {v\ cos9'\/m)~^ \ (d9/dy')y |, 
6 and 6' are the incoming and outgoing angles, and fi is the Maslov index given by the 
number of constant-energy conjugate points |Q. Performing the transverse integrations by 
stationary-phase approximation, valid in the h ^ limit applicable to the many-channel 



case that we are interested in, we obtain 55 



' ba 



V27iih 
2L,, 



sgn(a) sgn(6) 



' Ds exp { j^Ss{h, a, Ef) - i-fi. 



(C2) 



where the sum is taken over trajectories s between the cross section at x = with incoming 
and outgoing angles 6 and 6' such that sin 6' = dn/kLy and sin 6'' = bir/kLy (a = ±a, 
b = ±6). Since a and b are integers we see that the semiclassical approximation yields the 
intuitive result that only trajectories which enter and exit at discrete angles corresponding 
to the allowed quantized transverse momenta contribute to reflection. The reduced action 
(virial) is 



S{b,a,Ef) = S{yQ,yo,Ef) + fnvayo/Ly - h-Kbyy L 



y ' 



the new pre-exponential factor is 



mv I cos 6' 




(C3) 



(C4) 



and the Maslov index /2 is given by and the signs of the second derivatives of the action (see 
Ref. |^5[)- ^ similar expression holds for the transmission amplitudes, but the trajectories 
are now taken from the cross-section at x = to that aX x = L^. The Kronecker 5^ f, of the 
reflection amplitude, Eq. ( [7.3 b| ), has been taken care off by trajectories going from y to y' 
along the cross section x = (without entering the disordered region) not included in Eq. 
( |(J2D The rapid increase of the mean value of the reflection amplitude for high modes 
can be traced to the cos 6*' of Eq. (|C4|) since 6' approaches 7r/2 as a increases towards A^. 

When taking the average of the matrix elements of S over different impurity configura- 
tions, the long trajectories will have a rapidly varying phase which leads to zero average, 
while the short reflection trajectories which do not explore the bulk of the sample will sur- 
vive the average. In particular, trajectories corresponding to high mode numbers will have 
a large injection angle which makes very likely their reflection within a mean free path I 
from the entrance. The absence of self averaging to zero for the diagonal reflection elements 
subsists in the limit of L^; ^ oo of a semi- infinite strip, since it is a consequence of fast 
(short trajectory) processes. These short processes are the only ones contributing to {TrS) 
in the large A^-limit and this is the reason why in the numerics and in Eq. (|7.9|) the trace of 
the scattering matrix depends only on the properties at the entrance of the sample. How- 
ever, the reflected trajectories have a wide distribution of lengths, and we do not have as in 
Nuclear Physics ||5^ a clear-cut separation between direct (or fast) and long processes. 

In order to put these ideas in a more quantitative level, we calculate the semiclassical 
approximation to TrS"'. For n = 1 we only need the trace of the reflection matrices r and 
r', which can be evaluated with the aid of Poisson summation formula. 
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Tr r'''^' = ^raa = 7^ — E / '^^ *^^P C^^irnu) ^ sgn(M) sgn(u) 

a=l ^-^2/ m=-oo"'V2 ^(^g) 

VD^exp f -S'^(?2,M, E/) - , (C5) 

where u,u = ±u, not necessarily an integer. The stationary-phase condition for the integral 
over u is 

sgn(^) yQ - sgn(M) yo = 2mLy . (C6) 

Since < ?/o,l/o < Ly the only solution is m = m, ?/q = yo, and m = 0, i.e., the trajectory 
closes itself at the entrance cross-section with a specular reflection, where the transverse 
momentum (in the y direction) is conserved. Keeping only the m = term corresponds to 
the approximation of replacing the sum over a in Tr r^'^^ by an integral; the stationary-phase 
argument indicates that this is indeed the leading term in h. We then have 

Tr r-' = Y: exp fl.S,iyo,yo,Ef)-^^i^^ , (C7) 

V y s{yo,uo;yo,uo) 

Ss is the action integral along the classical path s, the stability prefactor is 



D. = D 



V nn^ I \ rlh^ I _ _ _ 

a=b=uo 



da? I , \ dl? 



(C8) 



and the new Maslov index fig is given by jls and the sign of the second derivative of the 
reduced action. For the semiclassical approximation to Tr r' we have a similar expression, 
given in terms of trajectories who have a specular reflection at the x = cross-section. 

Higher powers of S involve products of the matrices r, t, r', t' . For instance, the reflection 
submatrix of S"^ is rr + t't, and represents the two possible ways of returning to the entrance 
cross-section touching only once one of the vertical boundaries. The matrix products of 
the reflection and transmission amplitudes can also be obtained with the aid of Poisson 
summation formula, resulting in a stationary-phase condition very similar to Eq. ( |C6|) (but 
now the terms with yo and |/g come from actions of different trajectories, i.e. one from t 
and the other from if:') which indicates that we also have specular reflections at each bounce 
with the vertical walls. Once we take the trace we arrive at the intuitive result that TrS*" 
is given in a semiclassical approximation by closed trajectories, of a system with hard walls 
at the extreme cross-sections, with n bounces with the vertical walls. 

Bliimel and Smilansky |jl9| expressed TrS*" semiclassically in terms of periodic orbits of 



the Poincare scattering map, while in Ref. |^T[] it was shown that a periodic orbit of order n 
of the mapping corresponds to a periodic orbit of the closed system which hits n times the 
extreme cross-sections. This is precisely our statement of last paragraph. 

The semiclassical approximation to TrS**^, given by closed orbits hitting n times the 
vertical hard walls, is very useful since we can identify high values of n with long trajectories. 
The lower harmonics of p{6) are given by the ensemble average over short orbits, which is not 
zero in general; while the higher harmonics are given by long orbits leading to self averaging. 
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This explains why p{6) is relatively smooth and justifies the usefulness of our diagrammatic 
approach yielding only pi{6). 

According to Eq. ( |5.7] ) the two-level form factor is given by the ensemble average of 
iTrS""!^, therefore the semi classical approach also allows us to understand the phase shift 
correlations. The difference in the semiclassical approaches to the mean density and the 
TLFF is that the former is expressed as a one-phase sum while the latter a two-phase 
sum given by products of two periodic orbits of order n. In a diagonal approximation we 
would take both trajectories to be the same or symmetrically related. This yields a factor 
of two difference between the cases with and without magnetic field, which we obtain for 
small values of n (or in the slope at the origen) in our numerical simulations (see Fig. 2). 
Since the two curves for the TLFF must have the same area, they have to cross at one 
point, indicating that the diagonal approximation ceases to be valid for long trajectories. 
Therefore, for relatively low n we can approximate the Fourier components as simply 
given by periodic orbits hitting n times the vertical cross sections. 

The semiclassical approach allows us to give an intuitive interpretation of the Wigner 
time. According to Eqs. ( p.2| ) and ( |C2[ ) we have 



« J o-fi s{a,b) s'{a,b) 



(C9) 



The sum is over all pairs of trajectories s and s' contributing to r, r' , t or t' . In performing 
the ensemble average the short trajectories (where the diagonal approximation is valid) will 
dominate and for s = s' we have 

h ^ 

ME)) = -Y.T.Dsn, (CIO) 

ab s{a,b) 

giving the standard interpretation of the Wigner time as the mean scattering time (mean 
time of the classical scattering trajectories averaged with their stability prefactor) or the 
inverse escape rate |^ from the disordered region. 

The criterion developed Sec. ^ was to measure the departures from the circular en- 
sembles by the parameter t]c = TwEc/Tt, which is simply the ratio between the scattering 
and the diffusion times across the sample. The agreement with RMT obtained for large 
values of rjc translates into a Wigner time much larger than the diffusing time for traversing 
the structure (Ll/D). That is, the particle spends in the disordered region enough time to 
sample it before leaving. 

Separating the reflection and transmission trajectories s of Eq. ( P10|) we see that, at 



the semiclassical level, the Wigner time can be decomposed into transmission and reflection 
components. At the quantum mechanical level, this is not possible since the matrices r, r', t 
and t' are not unitary, and for a given impurity realization we obtain an imaginary part when 
using these matrices instead of 5* in Eq. ( |9.2D . However, these imaginary parts vanish under 
ensemble average giving well-deflned transmission and reflection Wigner times. These partial 
times can be related with the classical traversal times of a diffusion through a disordered 
region connected with leads (diffusion from an entrance to an exit p5| , ^ ). Calling Te = l/vf 
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the elastic scattering time, the quasi-one- dimensional case can be discretized as a diffusion 



in a one- dimensional lattice with m = L^/l sites. In this model it has been show |58.59 
that the mean traversal and reffection times are given by 

Tj, = Te m + ^) ' (Clla) 

l^ + l) ■ (Glib) 

Therefore, in the large m-limit, the mean traversal time is reduced by a factor of three with 

respect to the diffusion time in an infinite chain, ~ 1/3 Xg = 1/3 {L'^/2D). The 
classical transmission and reflection probabilities are 

T = , (C12a) 
m + 1 

R = — — , (C12b) 



and the total dwell time is 



Tr. 



TtT + TuR = Te m = Lr,/Vf , (C13) 



consistent with our findings that the Wigner time is approximately independent of disorder 
and coincides with the ballistic traversal time. For quasi-one-dimensional geometries the 
partial Wigner times obtained in the simulations after impurity average are in relatively 
good agreement with the average diffusion times r^T and r^i?. From Eqs. ( |(J11| ) and ( |U12[ ) 
we can also calculate the mean lenght of trajectories contributing to s^' , that is, those hitting 
n times the vertical walls at the entrance of the disordered region. For instance, for n = 2 
we have 

In the large m-limit, trajectories with n bounces will have a typical length proportional to 
n. 



APPENDIX D: INVARIANT MEASURE OF S IN THE POLAR 
DECOMPOSITION WITHOUT TIME-REVERSAL SYMMETRY 

In this appendix we find the metric of the scattering matrix in the polar decomposition 
for the case where there is no time-reversal symmetry. The calculation follows the same lines 



Al) we now have 



of that in Sec. |T|, the difference coming from the fact that instead of the decomposition 

S = X^S X = U' T U = Y' Y . (Dl) 
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The unitary matrix X diagonalizes 5", U and U' are independent block-diagonal unitary 
matrices given by the polar decomposition (|2.9|) . Y and Y' are independent unitary matrices 
that allow to define the metric of S. Diagonalizing F as in Eq. ( |3.4| ) we have the same Y as 
before (Eq. and 



u 



The infinitesimal variations of F, u*^^-*, and u^^-* are given by 

( |3.9| ) and ( |3.10D in terms of the antihermitian matrices 5Y ^ 5u^^^ and 5u^'^\ Analogously, 
we can write 

dY = Y' 6Y' , (D3) 

c^m(0 = y(0 5^(0 , 5^(0 = cia(0 + ^ cisW , / = 3, 4. (D4) 

Therefore the infinitesimal variations of 5* are given by 

dS = Y' {6Y' + 6Y)Y = Y' (idH) Y . (D5) 
The hermitian matrix dH defines the metric of 5* by 

2N 2N 

fi{dS) = n Re{dH)ij J] Im{dH)ij , (D6) 

i<j i<j 



and is given by 



dH-9 f ° VdQ-QdV\ 

[vdQ-QdV ] + 



+ 



In order to calculate the invariant measure of S it is convenient to break the product 



( p^ ) over the real and imaginary parts of the elements of dH as 



n Re{dH\a^, = n ^ ^l^^. dXa , (D8a) 



TV TV 



Y[ Im{dH)a,b Im{dH)a+N,b+N Re{dH)a,b+N 

(D8b) 
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Jl Re{dH)a,b Re{dH)a+N,b+N Im{dH)a,b+N Im{dH)h,a+N = Yl \ ? i " 1 — TT 11 ^-^2 ' 

a<b a<fe\^« + ^ ^& + ^//=l 

(D8c) 

n Re{dH),,a Re{dH)a+^,^^^ Im{dH),^,^, = [] {v', (^^'J + ^4'^) + Ql {ds^al + ^41^)) x 

a=l a=l 



The last product, Eq. ( |D8d| ), shows the redundancy of the polar decomposition in the 
absence of time-reversal symmetry. The matrix S is unchanged if in its polar decomposition 
we make the transformation [^| 



U^GU , U' ^ U'G^ , ^ = ( ^ ) ' 

where the blocks Q are pure-phase diagonal matrices, Qa = exp (irja). Under this transfor- 
mation we have 

ds^l^ds^l±dr^,, / = 1,2 (3,4) (DIO) 

and the r.h.s. of Eq. ( pSd] ) remains unchanged. Since drja are arbitrary we take them equal 
to ds^^^ 0] , and when multiplying both sides of the equation above by H^i ds^aa obtain 



the transformation of the volume element as 

n ReidHUa Re{dHU,,^^, Im{dHU^, ds^^} = J] 11 • (DH) 

a=l a=l V^a + -•■/ 1=1 

As in the time-symmetric case, we express the jacobian of the transformation in the Gibbs 
form ( [4.2| ), with an inverse effective temperature (3 = 2 and an effective hamiltonian ( [4.3|) 
given by the same logarithmic interaction as before (Eq. (f4.4|) ) and a one-body potential: 

y (A) = AT log (1 + A) , (D12) 

which coincides with ([4.5| ) up to order terms but differs in the terms of order unity (A^"). 
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TABLE I. Runs with zero magnetic field are indicated by R, runs with a nonzero flux 
through the sample are indicated by F (0o = hc/e is the elemental flux). NS is the number of 
samples (impurity configurations) considered. AR = L^/Ly is the aspect ratio of the samples and 
N is the number of propagating modes. Ef and k are the Fermi energy and the wavevector in 
Anderson units [Ef = —2 — 2 cos k). W \s the amplitude of the on-site disorder of the run. Tw is the 
Wigner time and A is the level spacing of a perfect (non-disordered) rectangle with sides Lx, Ly. 
The column showing -^-^^ ^ checks the approximate validity of Eq. ( |9.1D . I is the elastic- mean- free 
path and rc is the critical angle rfc (Eq. (|9.10D ) scaled with the mean phase shift spacing. (/>c is the 
slope of evaluated at rc obtained from Fig. 8. {g) is the averge dimensionless conductance and 
{Sg^) its variance. \(^rS)\ is the average trace of the scattering matrix. 
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FIGURES 



FIG. 1. Typical sample geometry where the disordered scattering is confined to the region 
< X < Lx of the wave guide. A and D are the amphtudes of the incoming flux, B and C are the 
outgoing amplitudes. 

FIG. 2. Phase shift form factor (Eq. ( [5.7D ) for quasi-ld conductors (upper inset) without 
magnetic field (Rl) and with (/>/(/'o = 2 (Fl). The large-M OF and UE values are given by thick 
solid and dashed lines. The number of impurity configurations is NS = 5000, the aspect ratio is 
AR = 4, the disorder is W = 1 and the number of propagating modes is = M/2 = 14. Lower 
inset: histograms of the phase shift distribution normalized by the uniform value po = M/27r. The 
Rl distribution is off-set by —1/3, while the Fl distribution is off-set by —2/3. 

FIG. 3. Phase shift form factor for quasi-ld insulators without (a) and with (b) mag- 
netic field. Filled circles correspond to strong localization, while diamonds represent a quasi-ld 
metal-insulator crossover (L^ ~ with weak disorder. Histograms in the insets show the corre- 
sponding phase shift distributions with an off-set of -1/3 for the strong disorder case and an off-set 
of -2/3 for weak disorder. The thick solid and dash lines in a (b) represent respectively the large-M 
form factor of one OF (UF) and of the superposition of two independent OF (UF). 

FIG. 4. Number variance, with (a) and without (b) time-reversal symmetry, as we approach 
the localized regime by increasing the disorder or the system length. Thick solid lines correspond 
to the number variance in the orthogonal (unitary) ensembles, Eq. ( |6.5| ), while thick dashed lines 
correspond to the number variance of the superposition ensembles, Eq. ( |6.7| ). White symbols show 
the approach to the quasi-ld localized regime by increasing the length keeping the disorder weak 
{W = 1): AR = 4 (squares), 10 (circles) and 30 (diamonds). Filled symbols show the approach 
to the localized regime by increasing the disorder keeping the geometry fixed {AR = 4), W = 2 
(triangles), 4 (circles), 6 (squares). 

FIG. 5. Diagrams of second-order perturbation theory for the average Green function used to 
calculate the average of the refiection matrix elements. 

FIG. 6. Average values of the diagonal reflection amplitudes as a function of the mode number 
for the samples scketched in the insets. The solid thick lines are guide-to-the eye with the functional 
form 1/ka and A^=14 modes, (a) Rl: AR = 4, W=l; R8: AR = 10, 14^=1; R2: AR = 4, W = 1.5. 
(b) R12: AR = 1, W=l, N=U; Rll: AR = 4, W=l; R13: AR = 1, W = 1.5, N = 14; R18: 
AR = 1, W=1.5, N=28. For the last sample we show the average diagonal reflection amplitude 
for only the first half of the mode numbers. 
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FIG. 7. (a) Phase shift form factors after unfolding the raw data for square samples (upper 
inset) for = 14 propagating modes. The sample with a disorder strenght of = 1 is represented 
by squares (R12) for B = and triangles (F12) for (j)/4>o = 2. Circles represent a sample with 
W = 1.5 and no magnetic field (R13). The large-M OE and UE values are given by solid and 
dashed lines. Lower inset: histograms of the phase shift distribution normalized by the uniform 
value po = M/2tt. The distribution R12 (F12) is off-set by -0.5 (-0.75), while the distribution 
of R13 is off-set by -0.25. (b) Phase shift form factors after unfolding the raw data (circles) and 
according to Eq. ( ^I7| ) (filled squares) for slab-shape samples (upper inset) without magnetic field 



(R19). Triangles: Form factors according to Eq. for (p/cpo = 2 (F19). The large-M OE and UE 
values are given by solid and dashed lines. Lower inset: histograms of the phase shift distribution 
normalized by the uniform value po = M/2tt. The R19 distribution is off-set by -0.4, while the F19 
distribution is off-set by -0.8. 



FIG. 8. Difference between the number variance obtained in the numerical simulations and the 
number variance of the standard ensembles. o"^(r) = $]^(r) — T,Q^{r) for samples without magnetic 
field and cr^(r) = S^(r) — S^^(r) for samples with magnetic field. Solid (dashed) lines are for 
samples without (with) magnetic field. Thin lines are for square geometries, R12, F12: AR = 1, 
W=l, N = 14; R13: AR=1,W = 1.5, = 14; R15: AR=1,W = 1.5, = 20. Thick lines are 
for thin slabs, R19, F19: AR = 1/4, W=l, N=U. 

FIG. 9. Slope of o"^(r) (obtained from Figs. 6 and 8) as a function of the inverse of the reduced 
critical angle Tc = {twEc/K) {M/2'k). The symbols used for the various runs are consistent with 
the labels of Table I of Appendix A. Samples Rl-Rll are quasi-one dimensional, while samples 
R12-R17 represent the two-dimensional case (square geometries). The dashed straight lines are 
guide-to-the-eye through the two groups of data points and illustrate the approximate validity of 
Eq. ( pup . 

FIG. 10. Energy dependence of the phase shifts for one of the quasi-ld conductors studied 
in Fig. 1 (aspect ratio AR = 4, disorder = 1) in an energy range where there are A'^ = 14 
propagating modes. 
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